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ABSTRACT 

We present the growing mode solutions of cosmological perturbations to the second order in the 
matter dominated era. We also present several gauge-invariant combinations of perturbation variables 
to the second order in most general fluid context. Based on the solutions we study the Newtonian 
correspondence of relativistic perturbations to the second order. In addition to the previously known 
exact relativistic/Newtonian correspondence of density and velocity perturbations to the second order 
in the comoving gauge, here we show that in the sub-horizon limit we have the correspondences for 
density, velocity and potential perturbations in the zero-shear gauge and in the uniform-expansion 
gauge to the second order. Density perturbation in the uniform curvature gauge also shows the corre- 
spondence to the second order in the sub-horizon scale. We also identify the relativistic gravitational 
potential which shows exact correspondence to the Newtonian one to the second order. 
Subject headings: cosmology: theory — large scale structure of universe 



1. INTRODUCTION 

Nonlinear evolution of perturbations is important to 
understand the observed large-scale structure in the uni- 
verse. Most of the study were made in Newtonian con- 
text (Peebles 1980; Vishniac 1983; Fry 1984; Goroff et al. 
1986; Makino et al. 1992; Fry 1994; Bernardeau et al. 
2002), but there have been growing interests in the cor- 
responding studies in Einstein's gravity context (Tomita 
1967, 1971, 1972; Kasai 1992, 1993; Bruni et al. 1997; 
Langlois & Vernizzi 2005; Malik & Matravers 2008; Ma- 
lik & Wands 2009; Langlois & Vernizzi 2010; Nakamura 
2010; Christopherson 2011). 

Our aim in this work is to present a complete set of 
growing mode solutions of zero-pressure scalar-type per- 
turbations to the second order. The solutions will be 
presented in a flat background with vanishing cosmolog- 
ical constant K = Q = K. Using the solutions we can 
clarify the relativistic/Newtonian correspondences of the 
density, velocity and potential perturbations available in 
various gauges to the second order perturbations. Previ- 
ously we have shown exact relativistic/Newtonian corre- 
spondence of the density and velocity perturbation vari- 
ables to the second order in perturbation (Noh & Hwang 
2004, NH2004 hereafter), which is possible only in the 
comoving gauge condition. In addition, the gauge trans- 
formation property of the metric and energy-momentum 
variables will be presented in the context of a general 
fluid with K and A. 

Section [5] is a summary of our convention and basic 
equations. In Section [3] we present the gauge transfor- 
mation properties and some useful gauge-invariant com- 
binations in general context. In Section |4] we show the 
growing mode solutions in various fundamental gauge 
conditions to the second order. In Section [5] we dis- 
cuss several distinguished relativistic/Newtonian corre- 
spondences. Section ini is a brief discussion. 

The reader who is interested only in our conclusion 
about the relativistic/Newtonian correspondences can go 



directly to Section [51 



2. SECOND-ORDER PERTURBATIONS 

We consider scalar-type perturbations to the second 
order in Robertson- Walker background. Our metric con- 
vention is (Bardeen 1988) 



ds^ = (1 + 2a) drf ~ 2a^/3,„dr;dx" 

{1 + 2^) gfp + J dx'^dx^, 



(1) 



where a is the cosmic scale factor, and a vertical bar in- 
dicates a covariant derivative based on . For the nor- 
malized fluid four- velocity we have Ua = —aw. a- To the 
second order, the basic perturbation equations of general 
fluicfl are presented in Equations (95)-(103) of Hwang & 
Noh (2007). In the case of a zero- pressure ideal fluid the 



^ For the energy-momentum tensor to the second order, see 
Equations (54), (84) and (85) in Hwang & Noh (2007). 
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equations become 



A 



K — 3Ha + 3(p-\ — ttX = "^0, 

. rr A + 3if 1 

47rGo/i + Hk ^ 05 = -rni, 



A + 3if 



X — 12nG{^ + p)av = 712, 



(2) 
(3) 
(4) 



k + 2Hk - AttG (S^ + 36p) + [3H+—]a = n3, (5) 



X + Hx — (fi — a — SttGH — 714, 

6il + 3H{6n + Sp) - {fi+p) (K-3Ha+ -Av 



= ns, 

[a\fi + p)vy 1 
a 



(6) 



(7) 



1 



a^{^i+p) a a{^+p) 

= "6, 



2A + 3if^ 



(8) 



where Ui are quadratic combinations of linear order per- 
turbations, and are presented in Equation (67)-(73) of 
Hwang & Noh (2007); k is the perturbed part of the 
trace of extrinsic curvature, x = a/3 + a^j is the shear 
of the normal frame vector, jj. is the energy density, p 
is the pressure, and H is the anisotropic stress; an over- 
dot denotes a time derivative based on t with dt = adi], 
and H = d/a. To the linear order, this set of equations 
is arranged in Bardeen (1988). The merits of this ar- 
rangement are (i) all perturbation variables are spatially 
gauge-invariant (this is true to the linear order, and can 
be extended to nonlinear order as well), and (ii) the tem- 
poral gauge (hypersurface or slicing condition) has not 
been taken, thus written in a sort of gauge-ready form. 
Setting any one of the perturbation variables equal to 
zero corresponds to a certain fundamental gauge condi- 
tion: see below. 

3. GAUGE TRANSFORMATION 

In this section we consider a general fluid in the pres- 
ence of K and A in the background. 

3.1. To the linear order 

We consider gauge transformation properties under 

x-=x'' + e(a;^) with e = e, e = and ^„ = i^/a; 
tildes indicate covariant quantities, whereas ^" is based 

(3) 

on 5„/3- 
To the linear order we have 

\a J a 
X^X-<\ K = K+(^3H+^^ae, (9) 

where a prime denotes time derivative based on rj. To 
the linear order, the following combinations are gauge 



invariant 

5y = S + 3{l + w) aHv = 3 (1 + w) aHvs, 
tfy = ip- aHv = -aHv^, tp^ = ip - Hx = -Hx^, 
11 

v^ = v x = — Xv, etc., (10) 

a a 

where w = p/fJ.. The gauge-invariant combination, for 
example, Sy is the same as S in the v = hypersur- 
face condition, thus 5y = S\y=Q. The temporal gauge 
condition, for example, v = fixes the temporal gauge 
mode completely. Thus, any perturbation variable in 
that gauge, for example S, can be equivalently regarded 
as a temporally gauge invariant ones, i.e., S\y = Sy. Sim- 
ilar complete gauge fixing is true for the following tem- 
poral gauge conditions: the comoving gauge (w = 0), the 
zero-shear gauge (x = 0), the uniform-curvature gauge 
(if = 0), the uniform-expansion gauge (k = 0), and the 
uniform-density gauge (^ = 0). The synchronous gauge 
(a = 0) is an exception, leaving a remnant temporal 
gauge mode C*(x) even after fixing the gauge condition. 

Similar complete gauge fixing and equivalence to 
unique gauge-invariant combination are possible for 
higher order perturbations as long as we take the spa- 
tial gauge condition in Equation ((T5|) : see Section VI of 
NH2004. Concrete construction of the gauge-invariant 
combinations to the second order will be presented be- 
low. 

3.2. To the second order 

The gauge transformation properties of the fluid vari- 
ables arcQ 



Sp = dp-p'^'^ -Sp'^'^ - 



-yee 



2 a 



a 

— V 

a 



H 



ai3 



H 



H' 



Q/3 



2^H 

a 



'2H^(q^ 



(11) 



(12) 



(13) 



(14) 



Gauge transformation properties of metric perturbations 
are presented in Equation (278) of NH2004. 

Now we address the issue of spatial gauge condition. 
The gauge conditions include the hypersurface (temporal 
gauge) condition and congruence (spatial gauge) condi- 
tion. To the linear order, the basic equations are already 
written in a spatially gauge-invariant form without di- 
rectly involving /3 or 7 separately (Bardeen 1988). All 
perturbation variables in Equations ©-(El) are spatially 

To the second order, the main fluid variables used in NH2004 
were based on the normal frame. The energy frame fluid variables 
are also present ed i n Eq uation (238) of NH2004 from which we 
have Equations I I11II - II14I I. 
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gauge-invariant or equivalently the same as the variables 
in the 7 = spatial gauge; this choice is unique because 
taking /3 = will leave a remnant gauge mode (Bardeen 
1988). Similarly, now to the second order we will take 

7 = 0, (15) 

as the spatial gauge (threading or congruence) condition. 
The choice is unique in the following sense. As 7 = 
completely fixes the spatial gauge mode, the remaining 



equations after fixing this condition can be equivalently 
regarded as spatially gauge-invariant ones even to the 
nonlinear order (Bardeen 1988; NH2004). Also our gauge 
related properties do not depend on our ignoring the vec- 
tor and tensor type perturbations which are coupled to 
the nonlinear order (NH2004). Under this spatial gauge 
condition we have ^" = even to the second order. 

From Equations pTjl - fM]) above, and Equation (278) 
in NH2004, we find the gauge transformation properties 
of the fluid and metric variables af0 



6p^ 6p p'e sp'e + p'ee + ^p"ee, 



I a 



Ila^ 


- 






a — 












V- 






a 


1 




+ 2 





a- 



1 /a" 



a a 
i— (^+ — 
a a 





1 /a" 


2 \~ 


1a- 
2 




1 

a 


{ax,c.e 



'2a + e)<\ 



-A- 



Ir'^e. " lex. - 3A-iv"v^ Qx,«e°, - lexp) 



t'2 



f 2a + ^ - 2e' - l-e) e.. + f « - 2^ - e' + -e- 

a \ 2 a J ' \ a a 



(16) 
(17) 

(18) 

(19) 

(20) 

(21) 

(22) 

(23) 



3.3. Gauge-invariant combinations 

From the gauge transformation properties of the metric and fluid variables we can construct the following gauge- 
invariant combination^ 



ify =ip — aHv — av 

1 /I 

+ 2 



^„ + 2Hip„ + -[H + 2H^] av 



^^ = ^-Hx-{^^ + 2H^^) X - i (ij + / + _L [^■'^^^^ _ A-i V" {xm)] 

3 1 A-iy^V' {xm) 



-i? A-iy" [2axx,a + (X - Hx) xA + -^a^HA'^ 



rX' X,a 



(24) 



(25) 



^ For K there is a typo in the sign of last term of Equation (279) 
in NH2004. 

For method, see Section VI. C. 2 of NH2004. The combination 
y)^ was presented in Equation (280) of NH2004; the other ones in 



NH2004 are based on the fluid quantities in the normal frame, thus 
here we are presenting again in the energy frame. 
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V>5 = V + 



+ 



+ 



3(1 + w) 
1 1 



1 . ^ ^ ^ H 1 il 



H 



2H'^ 2Hfi + w) 



-Xs + 



6{1 + w)H 



-Xs + 



3(l + «;) 
S 



6{1 + w)H 



(26) 



+ 



2a2 

A, 



2iJ/ 



1 /i 



/i 2 /i 



(5^ = 5 — — au — ( + ^—av 1 + — A ^V" (a^aw^a) , 



,Xv 



H 



2H J .a H 



6py = Sp — pav — ^Pv + 2^'^^ j + {(^v'^v^a) . 

^^=^~~^+~^~^'^" <^xX,a-x{aVyc)]„ + XX, 



rA- 



'X,a - 3A-iV"V^ (x 



1 . ^• 



,aX,/3)]| , 



-^A-1 
2 



A 



= K + ( 3-ff H — 2 ) ~ akyV + ^a^iiv^ — 2 (<^„ + aiJt;) — w + - (ay + ify + --aHv ) 



A 1 



a a 



-3aHA-'^ {ayAv + a^"v,a) ■ 
We introduced the entropic perturbation e as 

e = Sp— —S/j,, 

which is gauge-invariant only to the linear order. Wc can show that 



By = e — eav — \ p u 



■• ' ' 1 — av I av. 

fi 2 



(27) 
(28) 

(29) 

(30) 

(31) 

(32) 
(33) 

(34) 



Gauge-invariant combinations involving k are somewhat complicated. For example, we have 



.V A 



A 



3ir + 4 ) 



3H + 



A 



3H - 6HH - AH^ ^p^ 



-3 ( H^H + Ih'^ - ]-HH 



W^] +^(2^« + 3i^,.-2i^a.)^ 



1 



-—{H+-H 



rv Oi / • A TT TT \ — ^ 



Ak 



3H+^3H 



3H 



A\ 1 



-A-iV«VM-X«,ar^ 



1 K , 



2/ 2a2 



K,/3 



1 K'' 



Sir + 4 2 Sir + 4 3ij + 4 



3ij + ^ 2 3H+-^3H + 
I 



(35) 



where K/(3ir-|-A/a2) = [(3ii+A/a2)-iK;]. These gauge- 
invariant combinations will be used to relate solutions in 
different gauges in the next section. . 



4. GROWING MODE SOLUTIONS IN MATTER 
DOMINATED ERA 
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Now we consider a matter dominated era with K = 
= A. In this section we will present growing mode so- 
lutions of a zero-pressure fluid to the second order. Com- 
plete solutions to the linear order are presented in Hwang 
(1994). 



Thus ify — C (x) with vanishing transient solution where 
C (x) is an integration constant of the growing solution. 
To the linear order, the relatively growing solutions in 
expanding phase are 



4.1. Comoving gauge 
We take the comoving gauge 

w = = 7. 



(36) 



To the second-order perturbation, from Equations ([2]), 
(m and dH) we have 



'fv = ^ {no\v - n2\v) - ^"-el 
To the linear order we have simply 

ify = 0. 



(37) 
(38) 



(fv = C, ay ^ 0, Xv 
2 A 



2 1 
5H 



C, Ky 



2 A 



Sy 



c. 



(39) 



By plugging these linear solutions into rii of Equations 
©-([HI), we can obtain the full second order solutions: 
Equation ([5]) gives ay] Equations ^ and @ give ipy] 
Equation ^ gives Xv] Equation ^ gives Ky] and finally 
Equation ([3]) gives 5y. Then, we have 



1 



25 a^H^ 



^ ,a.i 



r 1 1 



ic'"C„ + A-V"(C„AC) 



= " + (^^^^ ~ 3A-iA-iV"V^ (CC^,)] - 



[2C'"C„ + 3A-1V" (C„AC)] 



i^(-2AC + 8CAC + 3C^"C„) ^ 
5 a^il 



1 



[2A(C'"C,c,) + 3V" (Cc,AC)], 



Sy = - 



175 a^H^ 

^ ^ (-2AC + 8CAC + 3C^"C,o)-h-4^^[A(C^"Ca) + 5V"(C,„AC)] 



5a2iJ2 



175 a4ff4 



(40) 



r 



These are a complete set of solutions in the comoving 
gauge. Using this set of solutions we can derive solutions 
in any of other gauge conditions. 



4.2. Zero-shear gauge 
The zero-shear gauge condition sets 



X = = 7. 



(41) 



This gauge is often referred to as the conformal Newto- 
nian or longitudinal gauge. To the linear order we have 



:C, a^ 



9 



A 



C. (42) 



Using these, we can calculate rii in the zero-shear gauge. 
Evaluating Equation in the comoving gauge we can 
derive ip^ to the second order. Starting from this. Equa- 
tion ^ gives a^; Equations ^ gives k^; Equation ([3]) 
and ([4]) give 5^^ and v^, respectively. Then, we have 



= + -C^ + — A-i [-CAC + 3A-iV"V'' {CC^a^p)] 



6 



1 



1 



25 
2 



175 a2i72 
6 1 
T75a2iJ2 



[C'"C,a -I- 5A-1V" (C,„AC)] 
[C'"C^a + 5A-1V" (C^cAC)] , 



,C + ;Jc2 + 4a-i [-CAC -f 3A-iV"V^ (^^C.„^)] 
5 25 25 



1 



175 a^H^ 
36 1 



[2C^"Ca + 3A-1V" (C,aAC)] , 



9 9 27 

^--5^^+50^^^+25^^"[-^^^ + ^^"^"^'(^^-'^)] ' 175a2i7 



[C^°C,„ + 5A-1V" (C,oAC)] 



6 



5v = 3 1 - ^ 



1 A 



1 1 



) C + + i^A-i [CAC - 3A- V"V'' (CC,,/^)] 



175 a^H^ 



[224CAC + 81C^°Ca - 36A~i V" (C ^AC)] 



1 



175 a4iJ4 



[A (C'"Ca) + 5V" (C„AC)]. 



(43) 



6 



HWANG ET AL. 



These are a complete set of solutions in the zero-shear 
gauge. Note that a^^ and Lp^ are presented in Boubekeur 
et al. (2009). 



4.3. Uniform- density gauge 
The uniform-density gauge takes 



7- 



(44) 



To the linear order, the solutions are 



Kg 



A 



15 a'^H'^ 



C, vs 



2 1 A 



1 

H 
2 1 



1 - 



15 a^H^ 



A 1 A 

-IjC, as — nrr9 fi 



XS = 7777 3 ^ 

15 H \ 



C. 



(45) 



To the second order, evaluating Equation ([26|) in the 
comoving gauge we can derive 



A 



15 a^H^ 



-CAC + -C'"C„ - A-^V" (C,aAC) 



75 o4 
1 1 



225 a^H^ 



^ 1^ [A (C'"C J + 5V" (C^AC)] + i (AC)' -f C^" (AC„) - A-^V^V^ [C„ (AC^)]| 
(AC'") ACc A-V^V^ [(ACa) AC^]} . 



(46) 



The rest of the solutions can be derived similarly: Equa- 
tion ([2]) gives Kg] Equation ([5]) gives a^; Equations ^ 
and ([7]) give X5 and vs, respectively. 



4.4. Uniform- expansion gauge 
The uniform-expansion gauge takes 

K = = 7. 



(47) 



It is often called as the uniform-Hubble gauge. To the 
linear order, the solutions are 



3iJ 



Xk 



3H 



14 2 1 3H 

-^C, <5k = -7 



3 A A 

^(j 



H 



5 3iJ 
3H- 



-C, 



2 1 1 A 



5aH3H 



rC, 



3 A 
5 ^ 



A 



rC. 



(48) 



To the second order, by evaluating Equation p5p in 
the comoving gauge or the zero-shear gauge we can de- 
rive </?K- The rest of solutions can be derived as follows: 
Equation ^ gives 6,^; Equation ([5]) gives a„; Equation 
([2]) gives Xk; Equation (|4]) gives As the solutions are 
lengthy we omit presenting them. Instead, we present 
some asymptotic behaviors. In the large-scale limit we 
have 



fK = C. 

In the small-scale limit, we have 



-a., = -a 



A 



A 



X' 

A 



Xk 



(49) 



(50) 



In the next section we will discuss these asymptotic so- 
lutions more. 



4.5. Uniform- curvature gauge 
The uniform-curvature gauge takes 

(p = = J. 



(51) 



It is often called as the flat gauge. To the linear order, 
the solutions are 



-C, Xv = 



c, K^ = m\-- + - 

3 1 



1 A 



C, w„ 



2 ■ 5a2i72 
1 



aR 



C. 



(52) 



By evaluating Equation (j27p in the comoving gauge or 
the zero-shear gauge we have 



2 A 



1 1 



4 1 



^C'"C„ - ^A-iy" (C„AC) 



175 a4ff4 



4CAC 

[A(C'"C„)-f 5V"(CaAC)]. (53) 



The rest of the solutions can be derived as follows: Equa- 
tion ([3]) gives Ktp', Equation ^ gives a^; Equation ([2]) 
gives Xf] Equation ^ gives vs- 

4.6. Large-scale evolution of curvature perturbation 

The curvature perturbations in various gauge condi- 
tions are known to be conserved in the large-scale limit. 
From our solutions, we have 

ffiv = fS = = C, 



12 
25~ 



25 



A-i [-CAC -f SA^iy^V'' (CC„/3)] ,(54) 



valid to the second order. General conservation proper- 
ties of ip^ , Lps and in the large-scale limit to nonlinear 
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order are studied in Lyth et al. (2005) and Hwang & Noh 
(2007). In addition, as we consider matter dominated era 
with fixed equation of state, Lpy^ is also conserved. Tfie 
following notations are often used in the literature 



H 



$yi and are often termed the Bardeen potential 
and Bardeen curvature, respectively (Bardeen 1980); 
these correspond to the Newtonian and the post- 
Newtonian gravitational potentials, respectively (Hwang 
et al. 2008). 

5. RELATIVISTIC/NEWTONIAN CORRESPONDENCES 

In this section we restore the speed of light c; in our 
convention the scale factor has a length dimension. In 
Newtonian context, the mass conservation, momentum 
conservation, and the Poisson's equation give (Peebles 
1980) 



1 



Vu 



a 

A 



1 



V • (5u) , 



1 



-u - Vu, 



$ ^ At:Gq5, 



(56) 
(57) 
(58) 



where (5, u, and $ are the relative mass density perturba- 
tion Sg/ g, the velocity perturbation, and the perturbed 
gravitational potential, respectively. By removing the 
gravitational potential Equations ([57]) and ([55]) give 

iv • (li + Hu) + AnGgS = -^V • (u • Vu) . (59) 
a 

In Newtonian context these equations are valid to fully 
nonlinear order. 



In Einstein's gravity, in the comoving gauge, from 
Equations (O and ([8]), Equations ^ and ([8|), and Equa- 
tion ([U, respectively, we have 



6y 



— S Y 



2Hkv - AnGgSy 



+ - 



Xv Xv,a0 



(Ax. 



C_ 

\2 



-Vu, 



(60) 

(61) 
(62) 

(63) 
(64) 



to the linear order. Equations (|60|) and (|61|) exactly coin- 
cide with Equations (f56| and (fSg]) . respectively. We have 
termed this coincidence, a relativistic/Newtonian corre- 
spondence of zero-pressure fluid in flat background to the 
second order in perturbation (NH2004, Noh & Hwang 
2005). Therefore, solutions for 5^ and valid for 

the Newtonian 5 and u, respectively, to the second or- 
der. 

From Equation (|58p with 5 = (5„ to the second order, 
we can construct relativistic combination of gravitational 
potential which exactly reproduces the Newtonian grav- 
itational potential to the second order. From Eqs. (jSj, 
(jlj and (lU, and using Equation (f28| . we have 



A 

K,v = ^Xv + n2\v 

By the identifications 

Sv = S, K, 

to the second order, and 
c. 



''^Xv = ^'^'^x = u = Vu, 



A-KGg ^ A A 

—dy + —ip^ = — 



A 



a' 

-1 



2ai 



2'KGg 



1 1 
— " ^x 



Yl-KGg 



X,a 



a -t- 



Ai^Gg 



(65) 
(66) 



where we have used Equations (|H), ([7]), ([5]) to the linear order. We have assumed only if = 0. Using these we can 
express the relativistic combination of the gravitational potential to the second order using Newtonian variables as 



2$2 _ 27rGga'^u^ - l27TGga^HA-^V ■ {u\76) - -A^^ (V$) 



c^a^ + \[- 2'KGga^u^ - UirGga^HA-'^^I ■ {u\75) - A" 
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(V$) + AnGga^u^ 



(67) 



Examination of the quadratic-order correction terms 
shows that those can be regarded as the post-Newtonian 
order corrections with ^ ~ K S order smaller 
than the leading order term, and become negligible in the 
small-scale limit. Therefore, with identifications oiSy, Ky 
and ^GR given in Equations (|63l) . (|64l) and (|67|) as the 
Newtonian density, velocity and gravitational potential 



fluctuations, respectively 

Sy=6, Ky = ~-\7-u, $GR = $, (68) 
a 

we have exact relativistic/Newtonian correspondence to 
the second order with Equations (|56p - ([58|) . 
For Newtonian gravitational potential, to the second 
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order, we have 



1 ^ 3 a^H^ , 



+ 



2 c2A 

_6 

175 a27J2 



20 



3 



(CAC) 



[C'"C,a + SA-^V" (C,aAC)] . (69) 



This, in general, differs from which also differs from 
—(fx as can be seen in Equation (1671) . However, to the 
linear order in general, or to the second order in the 



small-scale limit, we have 



-fx 



1$ 

2 ' 



(70) 



where a^^ and ip^ correspond to the Newtonian and 
the post-Newtonian gravitational potentials, respectively 
(Hwang et al. 2008). Therefore, in the sub-horizon limit, 
even to the second order we have 



Ki,(= Au^/a), 



-fx) 



(71) 



which correspond to the Newtonian perturbations (5, 
— iV • u and respectively. In the following we present 
further correspondences available in other gauge condi- 
tions, especially in the small-scale limit. 

In the small-scale limit, from our complete solutions 
we notice that 



A 

-fx 



Jip - 

A_ 
a 



2 c^A 
A 



C 



175 a^H^ 



2 c^A 



C 



tC 



6 c2 
' T75a2i72 



[A (C'"C„) -(- 5V" (C„AC)], 

[2A (C'"C,a) + 3V" (CaAC)] 



175 a^H^ 



[C'"Cc. + 5A-1V" (C,«AC)] 



(72) 



Thus, we have 



_ 2 c2A 
^ " 3a2ij2"^ 



(73) 



which is the relativistic version of Poisson's equation; 
it follows from Equations (g]) and ^ ioi K = 0. 
In Equation (|73|) instead of (5„ and ct y w e can replace 
other variable identified in Equation (|72p as well. We 
have shown that 6y and k„ exactly correspond to the 
Newtonian density and velocity perturbations, whereas 
= —fx corresponds to Newtonian potential pertur- 
bation in the small-scale limit. 

To the linear order, the perturbation variables S^, 
and axi= ~fx) ^''^^ known to correspond to the Newto- 
nian relative density perturbation, velocity perturbation 
and the gravitational potential perturbation, respectively 
(Bardeen 1980). Also to the linear order, all perturbation 



variables of density, velocity, and potential in the zero- 
shear gauge and the uniform-expansion gauge are known 
to have Newtonian correspondences in the sub-horizon 
limit: see Section 84 in Peebles (1980), and Table I in 
Hwang & Noh (1999). In this work we show that the 
same correspondences are now valid to the second order 
in the sub-horizon scale. We have the Newtonian corre- 
spondences for the following variable in the sub-horizon 
scale 



5 = 5, 
1 



A 



1 



A 

-V • u = — 

a a 

= = a. 



a 



-fx- 



(74) 



We can summarize the relativistic/Newtonian correspon- 
dences in a table form as the following: 





CG 


ZSG 


UEG 


UCG 


Exact: Second order 


Sy , K>y 








Linear order 


Sy , Ky 


"x(= -fx) 






Small-scale: Second order 


Sy , Hiy 


^x: Vx' "x(= -fx) 


5n, Vn, aK,{= -^K.) 


5ip 


Linear order 


Sy , K,y 


^x^ "x(== -fx) 


5k, Vk, aK(= -^k) 


5^ 



Table 1. Relativistic/Newtonian correspondences available in different gauge conditions. 

I 



The CG, ZSG, UEG, and UCG are acronyms of the 
comoving gauge, the zero-shear gauge, the uniform- 
expansion gauge, and the uniform-curvature gauge, re- 
spectively. Although the complete small-scale correspon- 



dences to the second order in the zero-shear gauge and 
the uniform-expansion gauge are impressive, lack of such 
correspondences in the other gauge indicates that the 
correspondence is still a non-trivial consequence avail- 
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able only under certain gauge conditions. 

We note that in our spatial gauge condition (7 = 0), 
the (temporal) comoving gauge {v = 0) does not imply 
the (temporal) synchronous gauge (a = 0). The subtle 
relations between our comoving gauge (u = = 7) and 
the original synchronous gauge (a = = /3) to the second 
order are studied in Hwang & Noh (2006). 

A complementary study of the correspondence between 
the linear perturbation theory and the post-Newtonian 
approximation is made in Noh & Hwang (2012). 

6. DISCUSSION 

In this work wc have presented the gauge transforma- 
tion properties and several widely used gauge-invariant 
combinations to the second order in the context of a gen- 
eral fluid: see Section [3] In the main part we have pre- 
sented the growing mode solutions of a zero-pressure fluid 
to the second order in several fundamental gauge condi- 
tions: see Sectional Based on the solutions we have clari- 
fied the relativistic/Newtonian correspondences available 
in different gauge conditions: see Section \E\ We have 
newly identified the gauge-invariant combinations of rel- 
ativistic perturbation variables which exactly reproduce 
the perturbed Newtonian gravitational potential to the 
second order: see Equation ([57| . Together with the 
previously identified gauge-invariant density and veloc- 
ity perturbations variables, now we complete the exact 
relativiatic/Newtonian correspondence to the second or- 
der perturbation: see Equation ([68]) . 

In Chisari & Zaldarriaga (2011) and Green & Wald 
(2011) the authors have considered situation where large- 
scale (near horizon-scale) fluctuations are linear whereas 
local density inhomogeneity is nonlinear. Both analy- 
ses were made in the zero-shear gauge (often termed the 
longitudinal gauge or the conformal Newtonian gauge). 
These two works are complementary to our present 
work in addressing relativistic identifications of New- 
tonian perturbations (the relativistic/Newtonian corre- 
spondences): their works cover large-scale linear per- 
turbation theory in the presence of small-scale nonlin- 
ear density inhomogeneity, whereas our work treats the 
second-order perturbations in all scales self-consistently. 
In the overlapping regime the results coincide. Both 
works have identified velocity and gravitational poten- 
tial perturbations in the zero-shear gauge {v-^ and a-^) 
and density perturbations in the comoving gauge (S^) 
as the corresponding Newtonian perturbation variables: 
this coincides with our identification to the linear order. 
For our extension of the identifications exactly valid to 
the second order, see Equation (l67|) . 

Extensions of the solutions in this work to cases with 
(i) general K and A, (ii) decaying solutions, (iii) gen- 



eral pressure, (iv) vector and tensor perturbations, (v) 
third-order perturbations, etc., are left for future studies. 
For our purpose of clarifying the relativistic/Newtonian 
correspondence our present study is enough (except for 
assuming A) as the correspondence in Equation (|67|) is 
valid up to the second order only ioT K = 0. Although 
the growing mode solutions in Section|3]are presented for 
A = 0, the correspondences in Equation (|67p are valid in 
the equation level with general A. 
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